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Ion mixtures require an exclusion core to avoid collapse. The Debye Hueckel theory, where ions 
are point charges, is accurate only in the limit of infinite dilution. The MSA is the embedding 
of hard cores into DH, and is valid for higher densities. The properties of any ionic mixture can 
be represented by the single screening parameter V which for the equal ionic size restricted 
model is obtained from the Debye parameter n. This Y representation (BIMSA) is also valid 
for complex / associating systems, such as the general n-polyelectrolytes. The BIMSA is the 
only theory that satisfies the infinite dilution limit of the DH theory for any chain length. 
Furthermore, the contact pair distribution function calculated from our theory agrees with 
the Monte Carlo of Bresme ea.(Phys. Rev. E 51 289 (1995)). 
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1. Introduction 

Stable ionic mixtures require an exclusion core to avoid the collapse of the system. 
For this reason the Debye Hueckel (DH) [ 1] theory, in which the ions are point 
particles, is accurate only in the limit of infinite dilution. The MSA [ El El E] that is 
an embedding of the DH theory into the hard-core Ornstein Zernike (OZ) equation 
is valid for high densities, and is even asymptotically exact at infinite density In 
the MSA the properties of any ionic mixture can be represented by a single screening 
parameter V [[6] which in the simplest equal ionic size restricted model is obtained 
from the Debye screening parameter k: 



r = — {vO + 2^-1} (l) 



k is the inverse of the Debye screening length defined by 
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where Zk is the electrovalence, pk is the density of component k, (3 is the Boltzmann 
thermal factor, cr^ is the hard core diameter, and e is the dielectric constant. The 
one parameter representation [ [7] is valid for a number of complex and associating 
systems such as dimers, and even polymers, where as in the Debye-Hueckel theory, 
the thermodynamic properties depend on a single screening parameter T: for the 
dimerization of equally charged hard ions we get 

<^=4r 2 (l + Ta) 2 (3 ) 
1 + 1 a 

After a detailed analysis[[H [9l [10] we can generalize eq.([3])in the form 

K 2 F(n,a) =4r 2 (l + To-) 2 (4) 

For dimers 

In the limit of total association a = 0, and we recover the DH limiting law with the 
charge of the polyelectrolyte, for any n > 3. Interestingly this will only happen if 
the correct virial (with the bridge diagrams) of the 3 body is used, and we get for 
linear homopolymer chains. Then after a lengthy calculation 

lim jF(n, 0) ~ n 2 ; correct DH limit 

p^O 



lim T{n, 0) ~ n; high density limit 

p^OO 



(6) 



The criticality of ionic systems, initiated by M.E. Fisher and collaborators is still a 
subject of current interest. [Qj], [12], [13]. In the original discussion of this problem 
Fisher and Levin [[11] used the DH [[I] theory in combination with Bjerrum associ- 
ation [I14j. It is clear that to get a mathematically well defined system all the ions 
need to have an exclusion core. Most of the subsequent work is directed at the in- 
clusion of the excluded core in the nonlinear Coulomb and the association problem. 
In our own previous work [ [T5] we used different combinations of the binding mean 
spherical approximation (BIMSA) [ [9J, [10] and various treatments of ion association, 
all derived from the work of Bjerrum [ 14J. The best agreement with computer sim- 
ulations was obtained from the SIS-BIMSA treatment of the association constant 
(Jiang, Blum, Bernard, Prausnitz and Sandler [OS]. In the SIS approximation of 
Stell et al [ [U] the chemical association constant is computed from the contact 
pair correlation function. A discussion of the merits of the different approximations 
was recently given by Aqua, Banerjee and Fisher [[12]). The real problem here is 
that, as it has been shown by Wertheim for the associating systems, the classical 
virial expansion does not converge [ [T71 [T8] and a new MSA (the BIMSA), based 
on the Wertheim-Ornstein-Zernike equation (WOZ) has to be formulated. As has 
been shown elsewhere [El ID], for dense systems, the scaling solution and the re- 
markably simple thermodynamics of the MSA applies to the new theories but with a 
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renormalized screening constant, T . In the ESMSA the low density limits are also 
included [ |21] . The real issue is that we need an internally consistent procedure and 
therefore there is an open (big) question about their reliability and accuracy that 
we will try to answer. 



2. Theory 

We wish a very happy birthday to Fumio Hirata on the occasion of his birthday 
Ever since the DH [H] theory was formulated, almost a century ago, there has been 
a steady effort in ways to improve and extend the range of its multiple applications. 
Most of the theoretical effort has been dedicated towards the inclusion of the hard 
exclusion core and the nonlinearity of the Poisson-Boltzmann equation. We know 
now that [ |5] that the MSA is asymptotically exact at very high densities. At lower 
densities the association limits yield exact conditions that are satisfied by closures 
of the Wertheim-Ornstein-Zernike equation: the exact DH limit at infinite dilution 
must contain the charges of the associated ions. The analytic solution of the MSA 
As has been pointed out in this work since the parameter space of these systems is 
large, we need to have an accurate, fully analytical theory for the properties of these 
systems. The general problem has been discussed in the past by several authors ( 
see for example [[19]). In recent work we have proposed a new theory, the ESMSA 
which could be such a theory, since there is only one screening parameter T, a in 
the MSA, but also the degree of dissociation, ionic diameters, effective dielectric 
constants enter in the calculation through the contact pair correlation function. 

The most interesting feature of the MSA even for complex systems, is that, 
as in the Debye-Hueckel theory, the thermodynamic properties depend on a single 
screening (scaling) parameter F: 

4[rf(l + ra) 2 = . 2 ^±^ (7) 

n is the inverse of the Debye screening length defined by 



An(3e 2 

k 



2 -*»h>c \ ^ 2 
« = — }^P kZ k 



and a is the degree of dissociation. The remarkable property of this equation is that 
it yields the correct asymptotic limits for zero density and also for very large density 
[E], where 

1 



- 2CT {Vr+2^-i>~|^-j ( 9 ) 

In the case of associating (polymerizing) ions the BIMSA for the restricted primitive 
model yields a very simple expression T B This result applies to Polyelectrolytes and 
explains the remarkable agreement with the simulations of Orkulas et al. [ [20] 
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The question is which is the best interpolation scheme between these two regimes? 
One possible answer can be obtained by testing different combinations as was done 
in our previous work [[15]: However it is clear that one cannot separate the sources 
of errors in the analysis of the criticality. However a more illuminating perspective 
can be obtained considering the simulations of Bresme et al [ 22J, as we will see 
it below. The usefulness of the MSA resides in the fact that the scaling solution is 
formally valid for the general mixture of arbitrary charge and size ions. 

It has been shown by Rosenfeld and Blum (Y. Rosenfeld and L. Blum,[|5] that the 
MSA is asymptotically exact at very large density. But the point remains that at low 
temperature ANY theory based on the normal Ornstein Zernike (OZ) equation will 
fail, because it will not be able to reach the full association limits. This is corrected 
by the Wertheim Ornstein Zernike (WOZ)[[T8] equation, which as has been shown 
to work remarkably well even for the limit of infinitely long polyelectrolyte chains [[8]. 

For many applications it is important to have an analytical, but at the same time 
an accurate theory of electrolytes. This is not an easy task, but recent advances have 
made this an attainable goal. One important ingredient of this theory is presumably 
the contact pair correlation function (PCF). In this communication we evaluate a 
simple functional form of the contact PCF, using the Monte Carlo simulations of 
Bresme et al. [ [22] . We find that the best analytical representation is obtained from 
this formula. It is increasingly erroneous at high charges. One possible reason for 
this failure is the fact that the simulations do not satisfy the PST, as is suggested 
by some preliminary simulation runs. Another source of errors is the neglect of the 
field parameter rj. 



3. The contact pair distribution function 

In the present contribution we show that the contact pair distribution function 
for a hard ions mixture can be represented explicitly by a functional of T [|9| [T0[, |2~T]: 

00/ \ HS, x vOl Wi _ *i 



(10) 

i] is a mean field parameter which depends on many body interactions, such as a, 
usually very small, and e is an effective dielectric constant. The hard core contact 
correlation is from the Carnahan-Starling approximation [ [2l] 

but in our case can be taken to be 1, since our system is very dilute. 

Indeed the proper discussion of this system is to treat it as a mixture of ions and 

dipoles for the associated part [[27]. We get for our restricted equal diameter problem 

&ij = <Ji = <Jj = a (12) 
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Figure 1. Comparison of the contact pair distribution function for 
gESMSA—BiMSA The squares are the Monte Carlo results of Bresme et al.[E2J. 
The triangles are the INV-C theory of Haymet et al The circles are the 

ESMSA-BIMSA with rj = and an effective dielectric constant of e = 1.17 
(which also depends on the association constant a) 



The simplest approximation is to take r] = 0, which means that we are neglecting 
higher order correlation effects [ [TO]. Then from eq. lfTO]) we get 

9™(<r) = e~'°V+r a) * (13) 

The results of the calculation are shown in figure 1. The agreement for ^- up to 
16 is quite good, and much better than the HNC theory [[20J which does not easily 
converge or the INV-C of Duh and Haymet [E5] which is shown in fig.l. However in 
our theory the error in the last point is sizable, and we have no explanation for this 
at the time, other than the suspicion that ths may be a consequence of the failure 
to satisfy the PST [|23]. 
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